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Abstract 

We show that a convex pyramid T in with apex at can be 
brought to the first quadrant by a finite sequence of monomial blowing- 
ups if and only if F fi (— IRq) — {^}- The proof is non-trivially derived 
from the theorem of Farkas-Minkowski. Then, we apply this theorem 
to show how the Newton diagrams of the roots of any Weierstrafi 
polynomial 

P(X, z) = + /li(x)z™-^ + • • • + hm~li^)z + h,ni^) , 

/ii(x) e /cjxi, . . . , a;ra][z], are contained in a pyramid of this type. Fi- 
nally, if n = 2, this fact is equivalent to the Jung-Abhyankar theorem. 

1 Introduction 

We will operate in the euclidean space with its affine structure. 
As it is classical, we will distinguish between the point-space X = R"' 
and the underlying vector space V = M". The vector addition is the 
canonical action (translations) of V on X. 

A polyhedron is the convex hull of a finite set of points generating 
the affine space X. Equivalently ([?], page 30), a polyhedron is the 
compact intersection of a finite set of half-spaces. Moreover, if ii^ = 
{Aq, Ai, . . . , Am} is a finite set of points generating X and, if 7^ = 
{Hi, . . . ,Hp} is the set of all (different) hyperplanes passing through 
all possible subsets of E consisting of n affinely independent points, 
then the set of vertices of the convex hull [E] of E is the set of points of 
intersection of all the subsets of TC consisting of n hyperplanes whose 
intersection is an only point^. 



"^cf., for instance, Vicente, J.L.: Notas sobre convexidad at http://www.us.es/da 
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A pyramid is the projection, from the origin G M", of a polyhe- 
dron contained in a hyperplane H not passing through 0. To be more 
precise: given a hyperplane H such that ^ H and a finite set E of 
points generating H, denoting by A the convex hull [E] of E, then 
the corresponding pyramid is 

r(A) = U (a)+ , 

aeA 

where (a)+ is the half-line of the non-negative multiples of a. Equiv- 
alently (c.f. Vicente, J.L., loc.cit.), a pyramid can be given by a poly- 
hedron A' in X, one of whose vertices is 0. In this case, the pyramid 
is nothing but 

r= U (a)+. 

a6A'\{0} 

Moreover, if {0,ai, . . . ,am} are the vertices of A', there is a hyper- 
plane H strictly separating from [ai , . . . , a^] ■ Then A = H f] A' is 
a polyhedron in H and T = r(A). 

Definition 1.1.— We will call a monomial blowing-up (resp. a mo- 
nomial blowing-down) any linear automorphism o/M" of the form 

(ai, . . . ,an) (ai, . . . ,an)Mij {resp. (oi, . . . ,a„) (ai, . . . ,an)Nij] 

where: 

1. i,j eZ, i^ j, 1 <i,j <n 

2. Mij (resp. Nij) is equal to the identity matrix in which the {i,j)- 
entry is set to 1 (resp. to —1). 

Remark 1.2.— With the notations of definition 11.11 the monomial 
blowing-up (resp. monomial blowing-down ) acts in the following way: 

(ai, . . . , On) (ai, . . . , + aj,...,an) 

resp. (ai, . . . , a„) (oi, . . . , Oj - Oj, . . . , a„) 

This corresponds to the behavior of the exponents of a monomial under 
the geometric monomial blowing-up Xi — > XiXj or the monomial blo- 
wing-down Xi — > Xi/xj. In fact, this geometric monomial blowing-up 
(resp. monomial blowing-down ) acts on monomials in the following 
way: 

•^1 ' ' ' "^j ' ' ' "^n -^l ' ' ' ■^i ' ' ' '^n 

ai an cti an 
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This is the reason of the name for these hnear automorphisms. 

From now on, we will use the name of monomial blowing-up (resp. 
monomial blowing-down) indistinctly for the linear automorphisms de- 
fined in or for the polynomial substitutions (resp.Xi — > 



Remark 1.3.— Let E = {ai, . . . , a^} be a finite set of points generat- 
ing a hyperplane H which does not contain the origin, let A = [E] C H 
be the corresponding polyhedron and r(A) the pyramid; then 



Equivalently, let A be the matrix whose row vectors are {ai, . . . , a„i}; 
then 



If (fij is the monomial blowing-up (resp. the monomial blowing-down) 
with matrix Mij (resp. Nij) and a- = ipij{ai) then E' = {a'^, . . . ,a^} 
generates a hyperplane H' not passing through 0. If A' = [E'], then 
r(A') = (p{T{A)), so it makes sense to speak on the transform of a 
pyramid by a monomial blowing-up or a monomial blowing-down. 

Definition 1.4.— The first quadrant of X is the seiMg. The opposite 
of the first quadrant of X is the set — Kq. 

The main problem we deal with in this paper is whether, given 
a pyramid r(A), it exists a finite sequence of monomial blowing-ups 
such that the transform of the pyramid by the sequence is contained 
in the first quadrant. We solve it by giving a geometrical criterion, 
from which we derive explicit computations using existing optimiza- 
tion algorithms. The criterion is the following: 

Theorem 1.5.— Let E = {ai, . . . ,am} be a finite set of points gen- 
erating a hyperplane H not containing the origin, let A = [E] be the 
corresponding polyhedron and r(A) the pyramid; then the following 
conditions are equivalent: 

1. There exists a finite sequence of monomial blowing-ups such that 
the transform of r(A) by the sequence is contained in the first 
quadrant. 



The second condition can be easily checked by the simplex method; 
in remark we will show how. Let us observe that the first condition 
implies the second because — Mg is stable by monomial blowing-ups. In 




r(A) = {(Ai, . . . , \m)A I (Ai, . . . , A„) G M^} . 



2. r(A) n {-w^) = {0}. 
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fact, if r(A) had a point a 7^ in common with — Mq then, no matter 
what sequence of monomial blowing-ups we apply, the transform of a 
will stay in — Kq. The point is then to prove that the second condition 
implies the first. 

As an application, we show that the Newton diagrams of all the 
roots of a WeierstraB polynomial 

P(x,z) = z"' + hi{x.)z"'-^ + - ■ • + /i„_i(x)2:+/i„(x) G fc[xl[z] , m > 1 

are contained in a pyramid satifying the equivalent conditions of the- 
orem [T31 Moreover, we show how, in dimension 2, this fact is in some 
sense equivalent to the Jung-Abhyankar theorem (c.f. [?]). 

2 The proof 

Remark 2.1.— Let us denote by A the matrix whose row vectors are 
{ai, . . . ,a.m}- We will speak of "bringing A to the first quadrant" as 
equivalent to bringing r(A) to the first quadrant. It is obvious that, 
if A has column with only positive entries, then it can be brought to 
the first quadrantby a finite sequence of monomial blowing-ups: it is 
enough to add this column to the others a suitable number of times. 

We use linear optimization methods to prove theorem ll.51 We refer 
to [?] for the theorem of Farkas-Minkowski and its consequences. In 
particular, we take from it (pages 42-51) the following consequence of 
this theorem (which might be also taken as an alternative statement 
of it): 

Corollary 2.2.— Let A be a matrix m x n; then one and only one of 
the following conditions hold: 

1. There exists a non-zero vector x > such that xA < 0. 

2. The system of inequalities Ay > has a non-negative solution. 
Inequalities must be understood componentwise. 

From this result we derive the following, which is the useful one: 

Corollary 2.3.— Let A be a matrix m x n; then one and only one of 
the following conditions holds: 

1. There exists a non-zero vector x > such that xA < 0. 

2. The system of inequalities Ay > has a positive integer solution 
(that is, a vector of positive integers). 
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Proof: By corollary 12.21 we must only prove that the existence of a 
non-negative real solution of Ay > implies that there is a positive 
integer one. Hence, it is enough to show the existence of a positive 
rational solution. 

Let y be a non-negative real column vector such that Ay > 0, 
let Oj be the rows of A, i = 1,... ,m, and let Uiy = Oi > 0. Let 
< e < minjaj} and let be a polidisc centered at y such that, for 
all z & K, one has \ — Oi + aiz\ < e. Then, —e < UiZ — Oi < e, so 
ai — e < aiZ < e + ai, which implies that aiZ > 0. Since K contains 
rational vectors greater than zero, the corollary is proven. q 

Proposition 2.4.— Let n > 1 and y = be a vector of 

positive integers with greatest common divisor equal to 1. There ex- 
ists a n X n matrix Y , whose entries are non-negative integers, with 
determinant equal to 1, one of whose columns is y. 

Proof: Let a G Z and denote by Eij{a), i ^ j, the n x n elemen- 
tary matrix, which is equal to the identity matrix with its (i, j)-entry 
replaced by a. Let us write y as a column vector. 

Let us assume that all the yj are multiple of one of them, say yi] 
then it must be = 1. Left multiplications by matrices Eji{—yj) allow 
us to transform the column vector y into a column vector having all 
entries equal to zero, except the i-th one which is equal to 1. Remark 
that all the elementary matrices we have used have a negative entry 
out of the main diagonal. 

Let us assume that no entry of the column vector y divides all the 
others and let be the smallest of all these entries. By assumption, 
there must be a j such that, in the euclidean division, yj = qjyi + rj 
with < r j < yi. Left multiplication of y by Eji{—qj) puts rj at the 
position j, leaving the other entries unchanged. In this way, we get a 
new column vector such that the greatest common divisor of its entries 
is 1 and the minimum of these entries has strictly decreased. Remark 
that, again, we have used an elementary matrix with a negative integer 
entry out of the main diagonal. If we repeat this process, after a finite 
number of steps, we fall in the preceeding situation. 

Summing-up: we have proven that, by left multiplication of the 
column vector y by elementary matrices having negative integer entries 
out of the main diagonal, we arrive at a matrix which is a column of 
the identity matrix. If we denote by Y the inverse of this product of 
elementary matrices, we see that y is a product of elementary matrices 
with positive integer entries out of the main diagonal, so 1" is a matrix 
with non- negative integer entries. Let In be the nxn identity matrix; 
then it is clear that Y = Yin contains a column equal to y and, of 
course, det(y) = 1. This proves the proposition. q 
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Proposition 2.5.— Let Y he a square matrix with non-negative inte- 
ger entries whose determinant is equal to 1. Then Y can he written as 
a product of (by order): a finite numher of monomial hlowing-up ma- 
trices, a permutation matrix and another finite number of monomial 
hlowing-up matrices. 

Proof: It is enough to show that Y can be brought to a permutation 
matrix by left and right multiphcaction by monomial blowing-down 
matrices. If q G then Eij{—q) = Eij{—lY; since Eij{—1) is a mo- 
nomial blowing-down matrix, then Eij{—q) is a product of monomial 
blowing-down matrices. 

Let yij be the smallest non-zero entry in Y . If all the elements 
in the j-ih. column are multiple of yij, we can get zeros in all the 
positions of this column, except (i, j), by left multiplication by ele- 
mentary matrices with a negative entry out of the main diagonal. In 
this case, the fact that det(y) = 1 implies yij = 1. If some entry in 
the j-ih. column is not a multiple of yij , say yij , and if yij = qyij + rij 
is the euclidean division, then left multiplication by Eji{—q) puts rij 
at the position {l,j) so the smallest non-zero entry of Y has strictly 
decreased. If we repeat this process, it is clear that we must arrive to 
the first case after a finite number of steps. The end of this process 
is a matrix Y( with a 1 at one position (denote it again by {i,j)) and 
zeros everywhere else in the j-th column. Moreover, Y( is the result 
of left multiplying Y by monomial blowing-down matrices. 

By symmetry, it is clear that we can get a new matrix Yi , obtained 
from Y( by right multiplication by monomial blowing-down matrices, 
and having 1 at the {i, j) position and zeros everywhere else in the i-th 
row and the j-th column. This is the basic argument of our proof. 

We may repeat the argument for the submatrix of Yi obtained by 
deleting the i-th row and the j-th. column, but seeing the operations 
in the whole Yi. This does not affect the form of Yi. The very end 
of the process is a matrix Yp, which is a permutation of the rows of 
the identity matrix, and which is obtained from Y by left and right 
multiplication by monomial blowing-down matrices. This proves the 
proposition. q 

Lemma 2.6.— In the situation of theorem \1.5l let A he the matrix 
whose row vectors are {ai, . . . , SLm}- For all x = {xi, . . . , Xm) € M™ \ 
{0} one has 



Proof: Let / : M" ^ M be a linear function such that /(aj) > 0, Vi = 
1, . . . ,m. There always exists such a function: we give an example. 



m 




i=l 
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Let g be an affine funcion such that H has the equation g = and 
g{0) < 0. Then, for ah i = 1, . . . , m one has = g{sLi) = g{0) + ^ {abi), 
so ~g{ai) > and we may thake / = ~g. 



which proves the lemma. □ 

Remark 2.7. Proof of theorem II.51 

Let us assume that is the only point of r(A) belonging to — Mq . 
By lemma for every x G \ {0}, the vector xA is different from 
zero. By assumption, it cannot be xA < 0. By corollary 12. 3[ there 
must exist a vector y € such that, written as a column vector. 
Ay > 0. We may assume that the greatest common divisor of the 
entries of y is 1. By proposition EIH there exists a matrix Y with non- 
negative integer entries and determinant equal to 1 such that y is one 
of its columns. Therefore, one of the columns of AY is Ay > 0, which 
implies by remark \2.1\ that AY can be brought to the first quadrant 
by a finite sequence of monomial blowing-ups. By proposition 12. 5( we 
obtain AY from A by applying to j4 a finite sequence of monomial blo- 
wing-ups, then a permutation of the columns and then another finite 
sequence of monomial blowing-ups. It is evident that the permutation 
of the columns plays no role: if the matrix with the permuted columns 
can be brought to the first quadrant, also the original one. This proves 
the theorem. 



The main application of theorem 11.51 we consider here lies is to the 
resolution of equations of the form: a WeierstraB polynomial equal 
to zero. Let k be an algebraically closed field of characteristic zero, 
X = (xi, . . . ,Xn) a collection of indeterminates, R = A;[x] the corre- 
sponding ring of power series, and let 



P(x,z) = z"' + hi{x)z"'-^ + --- + hm~i{x)z+h^{x) G klxj[z] , n>l 



be an irreducible WeierstraB polynomial; the object to study is the 
equation P{x,z) = 0. Let D E ^|x] be the discriminant of P with 
respect to z; the Jung-Abhyankar theorem (c.f. [?]) asserts that, if D 
is of the form x'*f7(x) with J7(x) G /s[x], U{0) ^ 0, then the roots 
of P(x, z) = are a full set of conjugate Puiseux power series in the 
variables x. Here, x*^ means x*^ = x"^ • ■ ■ x^", where a = (oi, . . . , a„) G 



Then, 




3 Applications 
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Zq \{0} When the discriminant has this very special form, we say that 
it is a normal crossing divisor. 

In general, the roots are not Puiseux power series in x. However, 
we can say something very important about them, namely 

Theorem 3.1.— The roots of P{x,z) = are power series belonging 
to a ring k{{x]-l^)) ■ ■ ■ {{xlfy)lx\^^ , . . . , a;!^**] whose Newton diagrams 
are contained in a pyramid r(A) such that r(A) fl (— Mq) = {0}. 

We will prove the theorem through several remarks. 

Remark 3.2.— Wc take the lexicographic order in the sense (c.f. [?], 
page 50): if a , b G then a <icx b if and only if the first component 
(from left to right) of a which is different of the corresponding in b is 
strictly smaller. 

Let us observe that a monomial blowing-up i^ij of the type a 
aMij with i < j preserves the lexicographic order, so it is an ordered 
automorphism of M** (endowed with the lexicographic order). In fact, 

^ijiflli ■ ■ ■ ) . . . , Oj, . . . , On) — (^1) • • • ) (^ii • • • , + Clj, • • • , On) ; 

so, if a <iex b, then: 

1. If <iex is decided before the position j, it is evident that (^(a) <iex 

2. If <icx is decided at the position j, this means that ai = bi, VZ, 
1 < ^ < J — 1 and aj < bj. Therefore, Oj + aj < bi + bj, hence 

99(a) <icx 99(b). 

3. If <iex is decided at a position I after j, this means that all the 
components of a until the (/ — l)-th coincide with the correspond- 
ing in b, so the same happens with <^(a) and (^(b). Since a/ < bi 
then 99(a) <iex (^(b). 

Wc call this an order-preserving monomial blowing-up. Notice that the 
corresponding monomial blowing-down ip^j^ is also order-preserving. 

Remarks 3.3.— Let A C Zq be a non-empty cloud of points; we 
call the transform of A by a monomial blowing-up (or a monomial 
blowing-down) the set of the transforms of all the points of A. 

3.3.1. Any monomial blowing-up 99 keeps Zq, that is, 99(Zq) c Zq. 
Therefore, if a , b G Zq are two points such that b G a -|- Zq then 
(^(b) G(^(a)+Z^. 

3.3.2. Let us assume that 99 = 99^^, p < q, is an order-preserving 
monomial blowing-up. Let a , b G Zq and write a' = 99(a), b' = 99(b). 
Let us assume that there exists an index j, 1 < j < n such that 
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bi > fli, Vi = then b'- > a[, Vi = In fact, the 

result is clear bv 13.3111 if j = n, so let us assume that j < n. Let c = 
(6i-ai, . . . ,6j-aj,0, ... ,0) and d = (0, . . . , 0, 6j+i-aj+i, . . . ,&„-a„); 
then b = a+c+d. Since a+c G a+Zg, then 93(a+c) e (/?(a)+Zo. On 
the other hand, since p < q, the first j components of f{d) are zero, 
so the conclusion is clear. It is obvious that the same happens if we 
replace ippq by the composition of a finite sequence of order-preserving 
monomial blowing-ups. 

3.3.3. For any j = 1, . . . , n and any v G Zq , we write Aj(v) = {v' G 
^0 \v'i> v^yi = then Ai(v) D ••• D A„(v). ByESEl 
for every order-preserving monomial blowing-up ip, every v G Zq and 
every j = 1, . . . ,n, one has (^(Aj(v)) C Aj((^(v)). Moreover, if $ is a 
composition of a finite number of order-preserving monomial blowing- 
ups, then $(Aj(v)) C Aj(^>(v)). 

3.3.4. Since 7^ A C , there is a minimum-lex u G A, so A C Ai(u). 
We will now prove that there exists a finite sequence of order-preser- 
ving monomial blowing-ups such that, calling $ the composition of all 
of them, one has $(A) C $(u) + . If A C An{u) there is nothing to 
prove, so we assume this is not the case. Let j be the smallest index 
such that A ^ Aj(u); then, necessarily j > 1. By the minimality 
of j, for every u' G A one must have u'^ > Ui, Mi = — 1 
Moreover, if u' G A \ Aj(u) then u'^ < Uj, so Uj > 0. Since u G A 
is the minimum-lex, for every u' G A \ Aj(u) there must exist an 
index i < j such that u'- > Ui. Let ii < j be the smallest index such 
that there exists u' G A \ Aj(u) satisfying > Ui^. Let <I>i be the 
composition of Uj monomial blowing-ups equal to 'Pi^j] then <I'i(u) = 
{ui , . . . ,Ui^, . . . ,UjUi-^ + Uj , . . . , Un) and, for every u' G A \ Aj (u) with 
u- > Ui, $i(u') = {u[, . . . ,u[^, . . . ,Uju'i^+u'j, . . . ,u'J. Since n^^ > Ui^, 
then Uj{u[_^ ~ ""ii) ^ '^j ^ % ~ u'j, hence Uju[_^ + u'j > UjUi-^ + Uj, so 
$i(u') G A,($i(u)). 

Now, ^>i(u) is the minimum-lex of <I>i(A) and, bv l3.3l3l $i(Aj(u)) = 
Ai($i(u)), Vi = 1, . . . ,n. ByESE ^i(A) C Ai($i(u)), Vi = 1, . . . , j- 
1. One must not forget that <I>i leaves invariant the first j — 1 compo- 
nents of every vector. Therefore, if <I'i(A) Aj($i(u)), there exists 
a smallest index 12 such that there exists u' G $i(A) \ Aj(<I>i(u)) 
satisfying u'^^ > Ui^. Necessarily 12 > h and we proceed as before, 
and so on. It is then clear that there exists a finite sequence of 
order-preserving monomial blowing-ups, whose composition $ is such 
that $(A) C Aj($(u)), Vi = 1, . . . , j. If there is still a ji such that 
<I>(A) ^ Ajj(<I>(u)), then ji > j and we proceed as before, and so on. 
This proves our assertion. 

Remark 3.4.- Proof of theorem 13.11 Let A be the Newton 
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diagram of the discriminant D of P{x,z); bv 13.3141 there exists a fi- 
nite sequence of order-preserving monomial blowing-ups such that, 
calling $ their composition, ^(A) C a -|- Zg where a E $(A). We 
make these monomial blowing-ups to act upon P(x, z) and denote 
by Q{x,z) the transform of P{x,z) by The discriminant D' of 
Q{x,z) is just the transform of D because D is a polynomial in the 
coefficients of the equation. Moreover, D' is a normal crossing divi- 
sor, hence the roots of Q = are all ordinary Puiseux power series, 
say with common denominator p of the exponents, because every irre- 
ducible factor of (5(x, z) has a discriminant which is a normal crossing 
divisor. If we come back to the beginning by applying the correspond- 
ing sequence of monomial blowing-downs, the region containing the 
Newton diagram of the roots of Q = 0; namely the first quadrant, 
obviously goes to a pyramid r(A) such that r(A) n (— Mq) = {0}. 
Since all the monomial blowing-ups are of the form ipij with I < j, 
we denote by i the minimum of all the indices I of these monomial 
blowing-ups, then $ leaves invariant the first i coordinates of ev- 
ery point, so the same happens with Therefore, every mono- 
mial x*^^^^ ■ ■ ■ x'^''^^x'^^-^^^ ■ ■ ■ x^^^ occurring in a root of Q evolves 
in a way such that the exponents ai/p, . . . ,ai/p remain unchanged. 
Therefore, if we fix a root of Q = 0, fix ai/p, . . . ,ai/p and write 
g' = x^i/P . . . xfPg"ixT+\'^^, xl-l^), with (J'^x'^C^l'^ . . . , xl-l^) € 
^I^i4a! ■ ■ ■ 1 ^"''^1) fo'^ the sum of all the terms of the root whose mono- 
mials start by x"^^"^ ■ ■ ■ a;"'"'^, the transform of £»' by produces a 

power series x"^^^ • • • a;^"'^^)']' ■ ■ ■ i ^iT^^) where q'[ has possibly 
negative exponents. Since • Zq is lexicographically well-ordered, 
so it is • Zq), hence the transform of q by belongs to 

k{{x\!^)) ■ ■ ■ {{x\^i))[x\^^ , . . . , which proves the theorem. 

When n = 2 there is much more to say, namely: 

Remark 3.5.— In our joint paper (cf. [?]), we prove the following 
for n = 2: 

1. The theorem 13. II without using the Jung-Abhyankar theorem. 

2. The Jung-Abhyankar theorem from the fact that the Newton 
diagrams of the roots lie in a pyramid r(A) such that r(A) fl 
i-Rl) = {0}. 

This shows that the Jung-Abhyankar theorem in dimension 2 can be 
proven by linear algebra techniques, without having resource to more 
sophisticated algebraic material. Moreover, in this case, the Jung- 
Abhyankar theorem is equivalent to the fact that the roots of the 
equation lie in a pyramid satifying the conditions of theorem 11.51 
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4 Short remarks on computations 



The explicit computations are a consequence, more or less obvious, 
of the convex calculus and the optimization of a linear function on a 
polyhedron by the simplex method. 

The point of departure will be always the list of points {ai, . . . , a^}, 
all different from 0, generating, either a hyperplane not passing through 
the origin, ot X = M". We add to the list, and write E = {0,a.i, . . . , am}; 
in both cases E generates the whole afHne space. We denote by A the 
matrix whose row vectors are {ai, . . . , a^}. 

Remark 4.1.— By elementary linear calculus (c.f. Vicente, J.L., loc. 
cit.), the (ra — l)-dimensional faces of the polyhedron [E] are produced 
by the following algorithm: we pick all the subsets of E consisting 
of n affinely independent points, and find the hyperplane determined 
by them; then we drop repetitions and keep only those hyperplanes 
leaving all the points of E in an only half-space. This algorithm is not 
the best possible, but improvements are out of the scope of this paper. 
If C = {-ffi, . . . , Hp} is the list of faces, then we get the vertices of 
[E] by the following algorithm: we pick all the subsets of C consisting 
of n hyperplanes whose intersection is an only point, find the point, 
drop repetitions and the remaining ones are the vertices. It is clear 
that E defines a pyramid T if and only if is a vertex of [E] . In this 
case, the faces of T are those Hi passing through 0. For instance, if 
we start from the points 

ai = (4, 7, -9) , a2 = (5, 7, -8) , ag = (3, 5, -9) , a4 = (4, 0, -1) , 

which generate Mq, the faces are 

7 xi — 13 X2 — 7 X3 = 18 — 9 X2 + X3 = 

-7x1 - 27x2 - 28x3 = 5x1 + 33x2 + 20x3 = 
19-2xi+X2 + 2x3 = 32 - 7x1 + 5x2 + 4x3 = 

and the vertices are E = {0, ai, a2, a3, a4}; therefore the points define 
a pyramid. The faces are normalized in the sense that all the points 
of E make their linear equations > 0. The faces of the pyramid are 
the first four and the edges are the positive half-lines determined by 
the four given points. 

Remark 4.2.— It is not difficult to know whether the pyramid F 
satisfies the condition F n (— Mq) = {0} or not. Let A = (Ai, . . . , A^) 
be a row of variables and let nii be the element in the i-th. column of 
the matrix AA; then F n (— Mq ) ^ {0} if and only if there is a feasible 
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solution to the set of linear constraints 

m 

l + ^^mi = 0, mi<0, Aj > , i = l,...,m. 

1=1 

The existence of a feasible solution can obviously be decided by the 
simplex method. In the preceeding example, the feasible solution does 
not exist, so F n (-M^,') = {0}. 

Remark 4.3.— It is also easy to find a positive solution of the system 
of inequalities Ay > 0, where ?/ is a column of variables. If rrii is the 
element in the i-th row of Ay, we can easily get a positive solution of 
Ay > by minimizing any of the coordinate functions on the set of 
constraints nii > 1 ,yi > 1, i = 1, . . . ,m. In our example, minimizing 
yi by the simplex method will produce the point (1,7/5,1), so the 
integer solution (5,7,5). 

The remaining computations to bring A to the first quadrant are 
straightforward matrix operations. 
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